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Abstract 



In order to obtain quite precise information about the shape of the particle paths 
below small-amplitude gravity waves travelling on irrotational deep water, analytic 
solutions of the nonlinear differential equation system describing the particle motion 
are provided. All these solutions are not closed curves. Some particle trajectories are 
peakon-like, others can be expressed with the aid of the Jacobi elliptic functions or 
with the aid of the hyperelliptic functions. Remarks on the stagnation points of the 
small-amplitude irrotational deep-water waves are also made. 



1 Introduction 

The purpose of the present work is to investigate the fluid dynamics at the 
propagation of small-amplitude gravity waves over irrotational deep water. It 
was widely believed that as small-amplitude gravity waves propagate on the 
surface, the particles of the fluid move on circular orbits, the diameter of which 
decreases with depth (see, for example, [TB], [30], [35], [31]). Indeed, in the 
first approximation, the fluid particles follow closed paths, but, by a phase- 
plane analysis of the nonlinear system describing the particle motion, in |12) 
and [8] it is shown that no particle trajectory is actually closed, unless the free 
surface is flat, there exists a forward drift over a period which decreases with 
greater depth. Similar conclusions hold also within the full nonlinear framework 
of the water-wave problem for symmetric periodic steady gravity waves (Stokes 
waves) travelling over a flat bed or over water of infinite depth. For a qualitative 
description of the particle trajectories in this framework, by methods from the 
theory of harmonic functions, see [6] and [21] (see also the very recent book 
[7]). The results in [6] are recovered by a simpler approach in [14j . where 
there are also described all possible particle trajectories beneath a Stokes wave. 
The particle trajectories change considerably according to whether the Stokes 
waves enter a still region of water or whether they interact with a favorable or 
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adverse uniform current. Some particle trajectories are closed orbits, some are 
undulating paths and most are looping orbits that drift either to the right or to 
the left, depending on the underlying current. 

Notice that there are only a few explicit solutions to the full nonlinear water- 
wave problems. For periodic gravity water waves in water of infinite depth, 
Gerstner constructed an explicit solution^ in 1802 [T5]. Gerstner's wave is a 
two-dimensional wave given in the Lagrangian description, by following the 
evolution of individual water particles. The motion of the water body induced 
by the passage of Gerstner's wave is rotational, it occurs in a flow with a specific 
non-constant vorticity. The fact that this flow is very special is confirmed also 
by the fact that this is the only steady flow satisfying the constraint of constant 
pressure along the streamlines cf. [31]. Beneath Gerstner's wave it is possible 
to have a motion of the fluid where all particles describe circles with a depth- 
dependent radius @], [25] . 

In order to obtain quite precise information about the shape of the parti- 
cle paths below small-amplitude gravity waves travelling on irrotational deep 
water, in this paper we provide analytic solutions of the nonlinear differential 
equation system describing the particle motion. We show that all these solutions 
are not closed curves. In the study of the nonlinear system ([33)) . a peakon-like 
trajectory (|54l) comes up. This solution has a vertical asymptote in the positive 
direction (see Figure 2). A peakon-like solution appeared also as particle path 
below small-amplitude periodic gravity waves travelling on a constant vorticity 
current (see [28] and [29] ) . The other solutions of the nonlinear system (|35|) are 
given by (|69l) . They cannot be expressed in terms of elementary functions. In 
some cases these solutions can be expressed with the aid of the Jacobi elliptic 
functions, in other cases they can be expressed with the aid of the hyperellip- 
tic functions. We draw some of the curves obtained for different values of the 
parameters (see Figure 4, Figure 5, Figure 6). We observe that some solutions 
(see ((89]) ) have vertical asymptotes in the positive direction (see Figure 6 too). 
This surprising feature is also found at the solutions describing the motion of 
the particles beneath small-amplitude capillary-gravity waves which propagate 
on the surface of an irrotational water flow with a flat bottom (see [27]). The 
particle seems to be shot out from the flow, this feature could reflect the wave- 
breaking phenomenon. 

For additional information on the particle paths within different types of pro- 
gressive water waves, in the framework of linear theory or in the framework of 
full nonlinear theory of periodic symmetric waves, and in the presence or not 
of the background currents and vorticity, see the following references: |5]-|14j. 
m~ EH, M- E2], [36], 04]. 

Another remarkable feature that we investigate in this paper are the stag- 
nation points for our problem, that is, points where the vertical component of 
the fluid velocity field is zero while the horizontal component equals the speed 
of the wave profile. The stagnation points are of special interest because they 



1 This solution was independently re-discovered later by Rankine |39) . Modern detailed 
descriptions of this wave are given in the recent papers [4] and 1231 . 
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are points where the flow characteristic often change. They could be located 
on the free surface, in this case the wave is called extreme wave, on the bot- 
tom or inside the fluid domain. For irrotational flows, the existence of extreme 
waves was predicted by Stokes [H], who also conjectured in 1880 that their 
profiles necessarily have corners with an angle of 120° at the crest. A rigorous 
proof of this conjecture has not been given until 1982, when it was established 
independently in [Tj and [38]. For rotational flows, the existence theory of small- 
amplitude waves is due to [T7] , [3D] , [35] , while the existence of rotational waves 
of large amplitude (approaching flows with stagnation points) in flows without 
stagnation points was recently established in [T3], [IS]. Some considerations 
about the possible shape of the limiting flow of the waves that almost admit 
a stagnation point are made in the paper [43] ■ But in the rotational case, for 
small-amplitude water waves with constant vorticity the stagnation points can 
also occur inside the fluid domain, as shown in [TB] and [35]. The papers [H] 
and [IT] contain existence results for small-amplitude steady waves with con- 
stant vorticity in the presence of stagnation points in the flow. The stagnation 
points around which the fluid rotates are centers. In irrotational deep water, 
there were revealed in [35] , although numerical and not completely rigorous (see 
[5J), stagnation points inside the fluid domain with discontinuous streamlines 
near the wave crests. These stagnation points are of saddle-point type. The 
formation of the stagnation points inside the flow may also be connected with 
a wave-breaking phenomenon in deep water; the traditional criterion for wave 
breaking is that horizontal water velocities in the crest must exceed the speed 
of the crest, see [2] where there are also reviewed different criterions, methods, 
measurements and laboratory studies to detect breaking events in deep water. 
Our solutions (|54|) and (|69| are also used to identify the stagnation points in 
small-amplitude irrotational deep-water waves. For the solution ([54"]) . a stag- 
nation point in the fluid appears only for t — ¥ ±oo. At this point the path of 
the particle has a horizontal tangent and the location of this stagnation point 
is on the bottom z = — oo. For the solution (|69[). the stagnation points are 
obtained by solving the equation (|TTj) . This equation can be solved graphically 
and depending on the signs and on the values of the involved parameters it can 
have one, two or three solutions. Which of these solutions are inside the fluid 
and their nature can be determined by a further study. 

2 The general problem for propagation of gravity waves in 
deep water 

Deep-water waves are modelled mathematically as periodic two-dimensional 
waves in water of infinite depth. To describe these waves we consider a cross sec- 
tion of the flow that is perpendicular to the crest line with Cartesian coordinates 
(x, z), the x-axis being in the direction of wave propagation and the z-axis point- 
ing vertically upwards. The water flow under consideration is bounded above 
by the free surface z — T](x,t) and has the bottom at z = — oo. The fluid is 
acted on only by the constant gravitational acceleration g, the effects of surface 
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tension being ignored. For gravity water waves, the appropriate equations of 
motion are Euler's equations (see |30j): 

u t + uu x + vu z = -- p Vx 
v t + uv x + vv z = -~p z ~ 9, 

where {u[x, z, i), v{x, z, t)) is the velocity field of the water, p{x 1 z, t) denotes the 
pressure, and p(x, z, t) is the density. Another realistic assumption for gravity 
water wave problem is the incomprcssibility, that is, the density p is constant 
(see [34]), which implies the equation of mass conservation: 

u x + v z = 0. (2) 

The boundary conditions that define water-wave problems come in various 
forms. For our problem, the dynamic boundary condition expresses the fact 
that on the free surface the pressure is equal to the constant atmospheric pres- 
sure denoted po- 

p = Po on z = r\(x,t). (3) 

The kinematic boundary condition expresses the fact that the same particles 
always form the free-water surface: 

V = r\ t + ur/ x on z = r)(x, t), (4) 

and the fact that at great depth there is practically no motion is expressed by 
the following boundary condition: 

(u, v) ->• (0, 0) as z -> -oo. (5) 

The governing equations and the boundary conditions ([T])-([5]) define the boundary- 
value problem for deep-water waves. 

An important category of flows is those of zero vorticity (irrotational flows), 
characterized by the additional equation: 

u z -v x = 0. (6) 

In what follows we will consider this type of flow. The idealization of irrotational 
flow is physically relevant in the absence of non-uniform currents in the water. 

For the water wave problem, there was a great interest in the regularity, or 
even stronger, the real-analyticity, of the solutions and particularly of the free 
surface. If these functions are real-analytic in some region, we can approximate 
them by using power series, approach used in small-amplitude wave theory. In 
the irrotational case, Lewy's theorem [33] shows that the free-surface must be 
a real-analytic curve and that the complex potential must have an analytic 
extension across the boundary (see [H]). Recent developments providing the 
regularity/real analyticity results for rotational flows were initiated in |10j for 
gravity water waves over a flat bed; for deep-water waves see [37]. In the case 
of zero vorticity, the approach in 10J yields an alternative short proof of the the 
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famous results obtained previously by Lewy. 




Figure 1 . Deep-water waves. 



3 Small-amplitude approximation of the water-wave problem 

Let us non-dimensionalise now the set of equations H])-©. If A > is the 
wavelength and a > is the wave amplitude (see Figure 1.), we make the 
following change of variables (for more details see [5U]). 



x H> Xx, z H> z, rj H> aij, 1 1-> -^=t 



(7) 



P i-> Po - gz + gp, (8) 

where, to avoid new notation, we have used the same symbols for the non- 
dimensional variables x, z, rj, t, u, v, p on the right-hand side. We set the 
constant water density p = 1. Therefore, in non-dimensional variables ©, ©, 
the water-wave problem ([I])-® becomes: 

u t + uu x + vu z = -p x 
pr(v t + uv x + vv z ) = -p z 
u x + v z = 

u z - -^v x = (9) 

v = a{rj t + ur) x ) on z = arj(x, t) 

p = ar/ on z — arj{x, t) 

(u, v) — > (0, 0) as z —} — oo. 

We observe in ((9]) that on z — ar) both v and p are proportional to a, this being 
consistent with the fact that as a — > we must have v — > and p — > (with no 
disturbance the free surface becomes a horizontal surface on which v = p = 0). 
For consistence, one requires that u is proportional to a too. Thus, by the 
following scaling of the non-dimensional variables: 



p H> ap, u i-> au, v H> av 1 



(10) 
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where we avoided again the introduction of a new notation, the problem ^ 
becomes: 

u t + a(uu x + vu z ) = -p x 
j?[v t + a(uv x + vv z )] = -p z 
u x + v z = 



II z 



X 



^v x = (11) 



v = rjt + aur] x on z = arj(x, t) 

p = i] on z = aij(x, t) 

(u, v) — > (0, 0) as z — > — oo. 

The classical approximation is the linearized problem obtained by requiring the 
amplitude of the free surface to be small, that is, by letting a — > 0, A being 
fixed. Thus, we obtain the linear system: 

u t + Px = 
j?vt + p z = 
u x + v z = 

hv x = (12) 



II T 



A 



V = rj t on z = 

p = X] on z = 

(u, v) — > (0, 0) as z-> — oo. 

From the third equation and the forth equation in (|T2"j) . we get that 



v zz + -^v xx = 0. (13) 



Applying the method of separation of variables, we seek the solution of the 
equation (|13|) in the form 



v(x,z,t) = F(x,t)G(z,t). (14) 
Substituting ([T4|) into the equation (p~3|) . we find 



+ < 15 > 

thus, 

G~dz^ __ A2F9?' ( ^ 

We observe in the above equation that the left hand side does not depend on z 
and the right hand side does not depend on x. Therefore, each side must be a 
constant, say 

F 9a; 2 G dz 2 A 2 1 J 

where K > is a constant that might depend on time. With the above choice, 
the solutions of the equations in (fl7|) are 

F(ar,t) = Asin(Kx) + Bcos(Kx) 
G(z,t) = Ce$ z +Ve-^ z , 
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where A, B, C, T> are constants depending on time. From the last condition in 
(TT2|) . that is, v — > as z — > — oo, we get that T> = 0, and thus, 



v(x,z,t) = Ce- z [Asm(Kx) + B cos(Kx)}. (18) 
On z = 0, by the fifth equation in (|12|). we have v = "q t , which yields 

C[^sin(iTa;) + S cos^a:)] = rj t . (19) 

Therefore, 

»(a?,z,t) = e^ z m. (20) 
Taking into account (|20p and the forth equation in (IT21 . we get 

z, t) = Y\ e ^ Zritx + J7<yX ' ^' ^ 21 ^ 

T{x,t) being an arbitrary function. The last condition in (|12p. that is, u — > 
as z — > — oo, yields 

T(x,t)=0. (22) 

The components u and u of the velocity have to fulfill also the third equation in 
(jT2J) . Replacing (|20p . (|2T|) with (|2^|) . in this equation, we obtain the equation 
that 77 has to fulfill: 

Vtxx + K 2 Vt = 0. (23) 
We seek periodic travelling wave solutions; thus, for the equation (|23|) with 

K := 2tt, (24) 

we choose the following solution 

r)(x, t) := cos(2tt(x - ct)), (25) 

where c represents the non-dimensional speed of propagation of the linear wave 
and is to be determined. With (|2"4l and (|25l) in view, the components of the 
velocity field become 

u(x, z, t) = ^^e^f z cos(27r(a; — ct)) 
A 

v(x, z, t) = 2nce^ z sin(27r(x - ct)). (26) 

In order to find the expressions of the pressure we take into account the first two 
equations in (fT2"j) and the expressions of the velocity field from above. Thus, we 
obtain 

2nc 2 

p(x, z, t) = ~^—e~ z cos(27r(x - ct)) + c (27) 
A 

where cq is a constant. On z = the pressure (|2"T|) has to fulfill the sixth 
equation of system (fT2j) . Hence, in view of (|25|) . we get 

2lTC 2 It, 

— — e~ z cos(2?r(a; - ct)) + c Q = cos(27r(a: - ct)). (28) 
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The above relation must hold for all values x £ R; therefore, we obtain 

co - 0, (29) 
and we provide the non-dimensional speed of the linear wave 

A 



c 



.2 



2vr' 



(30) 



We observe that the speed of the linear wave is proportional to the square root 
of the wavelength. In other words, long water waves travel faster than shorter 
waves. 

Summing up, the system fj 12[) has the solution 

rj{x,t) — cos(2tt(x — ct)) 
u(x, z, t) = ^]pe"? z cos(27r(x — ct)) 
v(x, z, t) = 2-Kce^ z sin(27r(a; — ct)) 
p(x, z, t) — e~r z cos(27r(a; — ct)), 

with c given by (l30l) . 

Taking into account 0, ©, (fTO)) . we return to the original physical variables. 
The speed of the wave (l30|) and the solution (|3T|) become: 

e - ± Ws- ± Vf (32) 

r)(x, t) = a cos(k(x — ct)) 
p(x, z, t) — po — gz + age kz cos(k(x — ct)) 
u(x, z, t) — acke kz cos(k(x — ct)) 
v(x, z, t) — acke kz sm(k(x — ct)), 



(33) 



k := f (34) 



where 



is the wave number. The sign minus in (|32[) indicates a left-going wave. We 
observe that the velocity components have the same amplitude acke kz which 
depends on position and decreases exponentially with the distance below the 
surface. 



4 Particle trajectories 



Let (x(t), z(t)) be the path of a particle in the fluid domain, with location 
(x(0), z(0)) :— (xq,zq) at time t = 0. The motion of the particles below the 
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small-amplitude deep-water waves with the velocity field (j3"3")l , is described by 
the following differential system: 

^| = u{x, z, t) = Ae kz cos{k(x - ct)) 

(35) 



dt 

where the constant A is 



^ - v{x, z, t) = Ae kz sin{k{x - ct)), 



A := ack ^ 0. (36) 



The sign of A depends on the sign of the wave speed c. Thus, if we choose in 
(|32[) the square root with minus, that is, we consider left-going waves, we have 
A < and if we choose in (|32|) the square root with plus, that is, we consider 
right-going waves, we get A > 0. 

To study the exact solution of the system (|3"5|) it is more convenient to re- 
write it in the following moving frame 

X = k(x- ct), Z = kz. (37) 

This transformation yields 



Integrating, we get 



If 



then, 



We denote by 



Z(t) = - log 



k A sin X(t)dt 



(38) 



^ = kAe z cos{X) - kc 
= kAe z sin(X). 
We write the second equation of this system in the form 

e - z dZ = kAsmX(t)dt. (39) 



er z = J kAsmX(t)dt. (40) 
J kAsinX{t)dt < 0, (41) 



(42) 



w = w(t):= J kAsinX(t)dt. (43) 
With (|41l) in view, we have 

to < 0. (44) 

From (|4"3")l we get 

kAsmX(t) = ^-. (45) 
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Differentiating with respect to t this relation, we obtain 

i a d 2 w 
kAcos(X)— = ^. 



dw 
~db 



From (|45[) we have furthermore: 

k 2 A 2 cos 2 (X) = k 2 A 2 



(46) 



(47) 



Thus, taking into account ([4"2"|). (ggj) , (gBJ and (g7]), the first equation of the 
system (1551) becomes 



a u) 



We make the following substitution 



A'r\/A' 2 .i- ( ^) • (48) 



(49) 



A being different from zero by (|36|) . Differentiating with respect to i this rela- 
tion, we get 

. d£ d 2 - 



^ dw dt 2 

We replace (J49J , ([50]) into the equation (|48l) and we obtain the equation 

aw if 

A solution of the equation (|5 1 [) is 

e = o, 

which, in view of (|4"!)| and f4"5)l implies 

sinX(t) = ±1. 



(50) 

(51) 

(52) 
(53) 



Therefore, from (|4"2")) with the condition (|4ip . and further from (j3"7| . a solution 
of the system S35\) is 



x(t) = ct + consti 

z(t) = -ilog(|fcAi+ consul) 



(54) 



const\ and consti are constants determined by the initial conditions (x(0), z(0)) 
(xo,^o)- We observe that 



lim = const, 



lim z(<) = lim z(t) = +oo, 



(55) 
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and 

lim x(t) — ±00, lim z(t) = —00. (56) 

Therefore, x = const will be a vertical asymptote and z = —00 will be a hor- 
izontal asymptote for the curve (I54[) . The particle seems to be shot out from 
the flow, this feature could reflect the wave-breaking phenomenon. For c > 0, 
the graph of the parametric curve (p>4"|) is drawn in the Figure 2. For c < 0, 
the arrows will change the direction. This peakon-like solution appeared also 
as particle path below small-amplitude periodic gravity waves travelling on a 
constant vorticity current (see [28] and [29]). 



Figure 2. Peakon-like trajectory. 



Calculating the derivatives of x(t) and z(t) with respect to t, we get 

x'(t) = c 



z '(t) = d 

V ' kA t-{-const-2 



(57) 



Hence, for the solution {5J$ a stagnation point in the fluid, where x'{t) = c, 
z'(t) = 0, appears only for t — > ±00. We observe that at this point the path of 
the particle has a horizontal tangent and is located on the bottom z = — 00. 
The other solutions of the equation (f5Tj) satisfy 

^ = 1 + kc. (58) 
aw w 



The homogeneous equation 



has the solution 



d 4 = (59) 
4 w 



£(w) = 9w, (60) 



where 9 is an integration constant. By the method of variation of constants, the 
general solution of the non-homogeneous equation (|58j) is given by 

£(w) = 6(w)w, (61) 



where 9{w) is a continuous function which satisfies the equation 
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With (|44l) in view, the solution of the equation (|62p is 

6{w) = fcclog(-w) + /3, 



(63) 



/3 being a constant. Therefore, the solution of the non-homogeneous equation 
(|58)l has the expression 



£(w) = w [fcclog(-uj) + 0\ . (64) 
Further, taking into account (l4l?|) . we get the equation that w(t) has to fulfill: 



^) = k 2 A 2 ~w 2 [kc\og{~w)+p] 2 



(65) 



From (|40|) and (|43|). between iu(i) and i?(i) there exists the following relation: 



w(t) 



-z(t) 



(66) 



Therefore, we obtain from (|65|) the differential equation that Z(t) has to satisfy: 



^) =fc 2 A 2 e 2Z - [/, rZ- >f 



By separating the variables, we get 



± 



dZ 



yJk 2 A 2 e 2Z -[kcZ - [3Y 



dt. 



(67) 



(68) 



Thus, taking into account (|39p and (f3"T|). i/ie oi/ier solution of the system H35\) 
has the following expression: 



x(t) = ct + jr arcsin 



A/1 



ct ± j: arcsin 



5 _Z(t)dZ(t)~ 




tit 





fccZ(t)-^ 



kAe z (f> 



(69) 



*(*) = i^W 

Z(t) being the solution of the equation A68\) . The constant /3 is determined by 
the initial conditions (x(0),z(0)) := (xo,z ). 

We observe that the solutions ()69[) are not closed curves. 
Indeed, if there exists ti > t\ such that Zfa) = Z(ti), then, we get zitz) — z(ti) 
and xfa) ~ x(t\) = c(t-2 — ti) ^ 0. If c > 0, the particles which follow these 
curves will have a forward drift, if c < 0, they will have a backward drift.D 



Let us make some remarks on the stagnation points inside the fluid. Calcu- 
lating the derivatives with respect to t of x(t) and z(t) from (|69[) . we get 
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With (|68p in view, for those Z satisfying the following equation 



kAe z 



kcZ- /3 



we have 



dZ d 2 Z 
~dt~' ~dP 
and thus, x'(t), z'(t) from (fTOf becomes 



At) 



z'(t) = 0. 



(71) 



(72) 



(73) 



Hence, for the solution &69\) , the stagnation points in the fluid are obtained by 
solving the equation |7_?| ). 

The equation (jTTj) can be solved graphically. Depending on the signs and 
on the values of the parameters k, c, A and j3, the equation (|71l) can have one, 
two or three solutions. See, for example, in Figure 3 some possibilities that can 



occur. With continuous line we have drawn 



kAe z 



Which of these solutions 



are inside the fluid and their nature can be determined by a further study. 



Figure 3. Graphical solutions of the equation (4.37). 

Let us now investigate more the equation (|68j) . Taking into account the 
expression of e 2Z as Taylor series: 

97 ^ (2Z) n „ 4Z 3 2Z 4 4Z 5 4Z 6 , . 

e^ = J2^-L = l + 2Z + 2Z^ — + — + — + — + ..• (74) 

n=0 y ' 



we get under the square root in (|68| the following power series 
(k 2 A 2 - p 2 ) + 2k(kA 2 + Pc)Z + k 2 (2A 2 - c 2 )Z 2 + 

4k 2 A 2 Z 3 2k 2 A 2 Z 4 4k 2 A 2 Z 5 4k 2 A 2 Z e 
+ 3 + 3 + 15 + 45 



(75) 



The constant A being different from zero (|36|) . the power series ([75]) contains 
for sure powers of Z higher than two. A partial sum of the above series is 
a polynomial of degree higher than two, therefore, in general, Z(t) cannot be 
expressed in terms of elementary functions. If in a partial sum of this series 
we consider powers of Z smaller or equal than four, then, the solution of the 
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equation (|68p involves an elliptic integral (for elliptic integrals see, for example, 
[3]). Its inversion would lead to an elliptic function. If in a partial sum of the 
series (|75|) we consider powers of Z higher than four, then, the solution of the 
equation (|68[) involves a hyperelliptic integral (for hyperelliptic integrals see, for 
example, [3 , page 252). Its inversion would lead to a hyperelliptic function. 

If in ([75)1 we consider powers of Z till three, then, the solution of the equation 
(|68|) involves the following elliptic integral of the first kind: 



±f-, ^ (76) 

j ^4k^v z3 + fc2 (2A 2 _ c 2)^2 + 2k(kA 2 + pc)Z + (k 2 A 2 - /3 2 ) 

This elliptic integral of the first kind may by reduced to the Legendre normal 
form. 

Case 1 : all the zeroes of the cubic polynomial under the square root in (|TG[) 
are real and distinct. We denote them by Z\ < Z2 < Z3. We introduce the 
variable tp by (see [40j Ch. VI, §4, page 602) 

Z = Z 2 sin 2 (f + Zi cos 2 tp, (77) 

and we get 

Ak 2 A 2 



(Z - Z X ){Z - Z 2 ){Z - Z 3 ) = 

sin 2 ip cos 2 ip(Z 2 - Z 1 f{Z 3 - Z x ) (l - k\ sin 2 tp) > 



3 

Ak 2 A 2 



3 

dZ = 2 sin tp cos ip(Z 2 — Zi)dtp, 
where the constant < k\ < 1 is given by 

« - l^f ■ ™ 

Therefore we obtain the Legendre normal form of the integral in (fTBf : 

1 f ^ =t, (79) 



Cl J ■ l\- k\ sin 2 if 



the constant factor in front of the integral being equal to 

1 



d :=±-j=k\A\ y /Z 3 -Z 1 . (80) 



The inverse of the integral in (|79|) is the Jacobi elliptic function sine amplitude 
sn (see, for example, [3]) 

sn (Ci t; k%) := sirup. (81) 
In view of the notation (f77|) . we get 

Z(t) = Z 2 sn 2 (Ci t;ki) + Zi cn 2 (Ci t;kx), (82) 
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cn being the Jacobi elliptic function cosine amplitude (see, for example, [3]). 
We introduce ([82]) in (|69|) and we get x(t) and z(t) explicitly. 

Case 2: the cubic polynomial under the square root in (1751) has only one 
real solution denoted Zq. We denote by p and q the real coefficients such that 

Ah 2 A 2 

^-Z 3 + k 2 (2A 2 - c 2 )Z 2 + 2k(kA 2 + pc)Z + (k 2 A 2 - p 2 ) = 
3 

Ah 2 A 2 

= ^-(Z-Z )(Z 2 +pZ + q). (83) 
We introduce the variable ip by (see 40 Ch. VI, §4, page 602) 



Z = Z + ^Z 2 + P Z + q tan 2 | , (84) 



and we get 

4fc 2 A 2 



., (Z-Z )(Z 2 +pZ + < ? ) 



4yfc 2 A 2 / /—„ \ 3 tan 2 f , , a 2 , 

\ / COS 2 



COS 2 



where the constant < k^ < 1 is given by 



■2 If, ^0 + 1 



2 I ^/Z 2 +pZo + q 



Therefore we obtain the Legendre normal form of the integral in ([76 

1 f dip 



C2j •/l-fc 2 sin 2 V 



t, (86) 



the constant factor in front of the integral being equal to 

2 



C 2 := ±^=k\A\(Z 2 +pZ + q)z. (87) 



The inverse of the integral in 

sn (C 2 t;k 2 ) = sinV> (88) 
Taking into account we get 
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We introduce ([59")) in ([fJ9")) and we get x(t) and explicitly. 

We observe that for that i's, denoted i + /C, with /C a period, for which the 

periodic Jacobi elliptic function cn satisfies the equation 

1 + cn (Ca t; k 2 ) = 0, (90) 

we have, in view of ([89]) and (|69|) . that 

lim :e(£) = const := a; + K , lim z(t) = oo (91) 
t-*(i+K) t-y{i+K) 

Therefore, x = x + K will be vertical asymptotes in the positive direction. This 
surprising feature is also found at the solutions describing the motion of the par- 
ticles beneath small-amplitude capillary-gravity waves which propagate on the 
surface of an irrotational water flow with a flat bottom (see [27]). The particle 
seems to be shot out from the flow, this feature could reflect the wave-breaking 
phenomenon. 

Let us draw below some of the curves obtained for different values of the 
parameters, using Mathematical. 

We consider k = 1, g = 9.8, a — 0.1 and (3 = 1. Then, by (]32]l . choosing 
the square root with sign plus, and by ([3"oT) . we get c = 3.1305, A = 0.31305. 
In this case, all the roots of the cubic polynomial under the square root in (|76[) 
are real and we get Z(t) in the form (|52"1) . The graph of the curve obtained is 
drawn in Figure 4. 




Figure 4. The solution for k=1, a=0.1, [3=1. 



For k = 2, g = 9.8, a = 0.1 and /3 = —1, choosing the square root with 
sign plus in ([52]). we get c = 2.21359, A = 0.442718, all the roots of the cubic 
polynomial under the square root in (f76j) are real, and thus, Z(t) has the form 
(|82|) . The graph of the curve obtained is presented in Figure 5. 




Figure 5. The solution for k=2, a=0.1, (3=-1. 



We choose now k = 4, g = 9.8, a = 0.1 and (3 — 1. We get by (l3"2"j) , choosing 
the square root with sign plus, c = 1.56525, A — 0.6261, but the cubic polyno- 
mial under the square root in (|76p has only one real root and we get Z{t) in the 
form ([59")) . The graph of the curve obtained is depicted in Figure 6. We observe 
that this solution has vertical asymptotes in the positive direction. 



2 In Mathematica the Jacobi elliptic functions are implemented as JacobiSN[«, m : = 
fc 2 ]:=sn(u; fci), JacobiCN[w, m := fc 2 ]:=cn(«; ki), JacobiDN[u, m := fc 2 ] :=dn(u; ki) 
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